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This lecture reviews current practice as well as new modeling ideas for the 
calculation of at least skin friction and heat transfer between the onset and 
end of transition. 

The principle means of transition zone modeling dating from the 
1950’s is through the intermittency models developed by Narasimha and 
colleagues. In these methods the length of the transition zone is closely 
correlated to mean boundary layer properties at the onset location and so 
the procedure is critically dependent on a correct determination of the 
onset location. The procedure is independent of the level or character of 
the free-stream disturbances implying independence of the spot 
generation mechanism or other detailed physics. Determination of onset 
location is principally by correlation methods of by e N for quiescent 
streams. Recent development of PSE methods by Herbert and others 
provides a powerful method for calculating onset at all turbulence levels. 
Additional data for modeling of bypass transition over heated flat plates 
have been obtained by Sohu and Reshotko in the NASA-Lewis boundary 
layer channel. With a series or grids placed upstream of the contraction, 
free stream turbulence levels ranging from 0.4% to 6% can be produced. 

For Tu=0.4%, heating the plate advances transition as would be expected 
fro T-S considerations. For Tu>l%, there is negligible effect of heating on 
transition. 

Rai and Moin (1990) have performed a direct Navier-Stokes (DNS) 
simulation of the grid 2 case of Sohu and Reshotko. This calculation which 
used about 8000 hrs. of CRAY-YMP time reproduced the onset of transition 
and most profile features through transition but missed the turbulent skin 
friction level. Such DNS simulation have to do better before they can be 
used as validation “experiment” for simpler modeling methods. 

Sohu’s intermittency conditioned mean profiles show that the 
nonturbulent part is Blasius near the beginning of transition but is fuller 
than Blasius as transition progresses. The turbulent part is initially less full 
than the standard turbulent boundary layer but approaches the Musker 
profile as transition is completed. Sohu’s boundary layer spectra for Grids 
2 and above are unambiguously turbulent. The spectra for Grid 0.5 and 
Grid 1 are somewhat ambiguous but our sense is that the bypass made is 
dominant. Sohu’s spectra have some t-s features. Suder’s spectra in the 
same facility for Grid 0 without plate heating fairly clearly document a T-S 
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path to transition. The turbulent spectra observed in bypass situations 
suggest that bypass transition might be modeled using turbulent flow 
methods. Additional experiments are under way for boundary layers with 
pressure gradients. 

The balance of the presentation shows the experience to date of the 
Wu and Reshotko with a multiple-time-scale k-e method. The skin friction 
results for the ERCOFTAC T3A case (Tu=2.8%) are reasonable but subject to 
improvement. The distributions of the k and e components through 
transition are shown. More specifically, the low frequency component of k 
(kp) through transition is comparable to what is expected from such 
experiments including law-of-the wall similarity in the near-wall region. 
For turbulence levels above 12%, the law-of-the-wall similarity breaks 
down indicating that the outside turbulence is strong enough to change the 
near wall behavior thus altering the characters of the turbulent boundary 
layer at sub-elevated turbulence levels. 

Work is continuing toward improving the MTS modeling and making 
the procedure more robust. Also it is important that the receptivity issue 
be examined toward identifying the physics of bypass initiation. Of 
particular interest in this regard are the forced algebraically growing 
subcritical modes identified recently by Henningson, Trefethen and 
colleagues. 
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Schematic dibgrak of the wind tunnel 
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Pig. 28(a) Interoittency profiles across the boundary layer for grid l 
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Fig. 29(b) Conditionally saapled sean velocity profilM in wA onlts 
at X-ll inches for grid 1 
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Fig. 29(c) Conditionally sampled naan velocity profiles in wall units 
at X-13 inches for grid 1 
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Pig. 29(e) Conditionally sampled mean velocity profiles in well units 
at X-17.5 inches for grid 1 
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Fig. 54 Freestteain spectra for 4 grids 
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Fig. 53(d) Boundary layer spectra measured at locations where rms velocity 
reached maximum for grid 2 
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Multi-Time-Scale Low Reynolds Number k-e Turbulence Model 
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Multi-Time-Scale Low Reynolds Number k-e Turbulence Model 
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INFLUENCE OF THE STARTING LOCATION 
Initial profiles: K and e are Rodi's profiles 
Freestream: K =0.7K„ 
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ERCOFTAC T3A 





ENERGY TRANSFER FUNCTION, 
Multi-time-scale k-£ Model 
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DISSIPATION, E T , 
Multi-time-scale k-s Model 
ERCOFTAC Test Case T3A 









TURBULENT KINETIC ENERGY, K , 

Multi-time-scale k-c Model 

High Free— Stream Turbulence, 15 percents 
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TURBULENT KINETIC ENERGY, 
Multi-time-scale k-£ Model 
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